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ABSTRACT 

en ! 

The identification of type Ib/c supernovae as GRB progenitors is motivated 
by the association of GRB980425 with SN1998bw and of GRB030329 with 
SN2003dh. While the 7-ray luminosity of GRB030329 was typical to cosmo- 
logical GRBs, the luminosity of the nearby (40 Mpc) GRB980425 was ~ 5 orders 
of magnitude lower. The large luminosity difference is commonly explained by 
hypothesizing that either SNe Ib/c produce two different classes of GRBs, or that 
I GRB980425 was a typical cosmological GRB jet viewed off-axis. In the latter sce- 

nario, strong radio emission, L v ~ 10 30 z/ 10 Q^ z erg/sHz, is expected at ~ 1 yr delay 

(N) \ due to jet deceleration to sub-relativistic speed, as observed from GRB970508. 

The radio luminosity of SN1998bw was 3 orders of magnitude lower than this 

t— 1 \ value. We show that the low radio flux may be consistent with the off-axis jet 

interpretation, if the density of the wind surrounding the progenitor is lower 
than typically expected, rh = (M/lO~ 5 M yr -1 )/(t> w /lO 3 kms -1 ) ~ 0.1 instead of 
P-f rh > 1. The lower value of rh is consistent with the observed radio emission from 

the supernova shock driven into the wind. This interpretation predicts transition 
to sub-relativistic expansion at ~ 10 yr delay, with current ~ 1 mJy 10 GHz 
flux and my « 23 optical flux, and with « 10 mas angular source size. It also 

.^J I implies that in order to search for the signature of off-axis GRBs associated with 

nearby Ib/c supernovae, follow up observations should be carried on a multi-yr 
& ; time scale. 

Subject headings: gamma rays: bursts and theory — supernovae: general — 
supernoave: individual (SN1998bw) — radio continuum: general 



1. Introduction 

The association of GRBs with type Ib/c supernovae is motivated by the temporal and an- 
gular coincidence of GRB980425 and SN1998bw (Galama et al. 1998b), and by the identifi- 
cation of a SN1998bw-like spectrum in the optical afterglow of GRB030329 (Stanek et al. 2003; 
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Hjorth et al. 2003) (See, however, Katz 1994, who suggests that a SN like emission may re- 
sult from the impact of the relativistic GRB debris on a nearby dense gas cloud). The 
7-ray luminosity L 1 ^ so ~ 10 51 erg/s of GRB030329, inferred from the redshift z=0.1685 of 
its host galaxy, lies within the range of typical cosmological GRBs, L 7) i so . — 10 52±1 erg/s 
(e.g. (Schmidt 2001)). The subscript "Iso." indicates luminosity derived assuming isotropic 
emission. The association of GRB980425 with SN1998bw sets the distance to this burst to 
38 Mpc (for H = 65km/sMpc), implying that its luminosity is nearly 5 orders of magnitude 
lower than that typical for cosmological GRBs (Pian et al. 2000). 

Two hypotheses are commonly discussed, that may account for the orders of mag- 
nitude difference in luminosity. First, it may be that SNe Ib/c produce two different 
classes of GRBs, with two different characteristic luminosities. It is now commonly be- 
lieved that long duration, T > 2 s, cosmological, L 7jIso . ~ 10 52 erg/s, GRBs are pro- 
duced by the collapse of SN Ib/c progenitor stars. It is assumed that the stellar core 
collapses to a black-hole, which accretes mass over a long period, ~ T, driving a rel- 
ativistic jet that penetrates the mantle/envelope and then produces the observed GRB 
(Woosley 1993; Paczyhski 1998; MacFadyen & Woosley 1999). This scenario is supported 
by the association of GRB030329 with SN2003dh, and by additional evidence for opti- 
cal supernovae emission in several GRB afterglows (Bloom 2003). The origin of a sec- 
ond, low-luminosity, class is unknown. It may be, e.g., due to supernova shock break-out 
(Colgate 1968; Woosley & Weaver 1986; Matzner & McKee 1999). The small radius and 
high density of carbon/helium SN Ib/c progenitors may allow acceleration of the shock to 
mildly relativistic speed as it propagates through the steep density gradient near the stellar 
surface. It is not clear, however, that the energy transferred to mildly relativistic material 
is sufficient to account for the 7-ray emission. 

A second possibility is that GRB980425 was a typical, cosmological GRB jet viewed 
off-axis (Nakamura 1998; Eichler & Levinson 1999; Woosley, Eastman & Schmidt 1999) 
(Granot et al. 2002; Yamazaki, Yonetoku & Nakamura 2003). Due to the relativistic expan- 
sion of jet plasma, with Lorentz factor T > 100 during 7-ray emission (Krolik & Pier 1991; 
Baring 1993), 7-rays are concentrated into a cone of opening angle ~ 1/r around the ex- 
pansion direction. Thus, if the jet is viewed from a direction making an angle larger than 
6*j+few xl/r with the jet axis, where 0j is the jet opening angle, the 7-ray flux may be 
strongly suppressed. In this scenario, strong radio emission, L v ~ 10 3< V 10 Q^ z erg/sHz, is 
expected at ~ 1 yr delay (Frail, Waxman & Kulkarni 2000; Livio & Waxman 2000) as the 
jet decelerates to sub-relativistic speed and its emission becomes nearly isotropic (Perna & 
Loeb (1998) have suggested that radio emission from the bow shock surrounding the jet 
may be detected on shorter time scale). Radio emission associated with transition to sub- 
relativistic expansion has been observed for GRB970508 (Frail, Waxman & Kulkarni 2000, 
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hereafter FWKOO), for which the transition was inferred to occur on ~ 100 day time scale. 
A flattening of the light curve, which is expected to accompany the transition (Waxman, 
Kulkarni & Frail 1998, FWKOO, Livio & Waxman 2000), has been observed on a similar time 
scale in the radio light curves of most well-sampled afterglows (Frail et al. 2003). 

Based on the radio observations of GRB970508 and on the hints for an association of 
GRBs with supernovae, Paczyhski (2001) suggested to search for radio emission from ~ 1 yr 
old "GRB remnants" among nearby (< 100 Mpc) supernovae. Levinson et al. (2002) have 
shown that a large number of such remnants may be identified by all sky radio surveys. A 
radio survey monitoring 33 type Ib/c supernovae for m 1 yr has recently been published by 
Berger et al. (2003b). The bright radio signature expected from (decelerated) GRB jets 
has not been detected, leading to the conclusion that the vast majority of type Ib/c SNe 
are not associated with cosmological (L 7) i so . — 10 52 erg/s) GRBs. The fact that the radio 
luminosity of SN1998bw after ~ 1 yr delay was 3 orders of magnitude lower than the expected 
L v ~ 10 30 z/ 10 QH Z er g/ s IPz, appears therefore to rule out the "off-axis jet" interpretation of 
the low luminosity of GRB980425. We show here that this is not necessarily the case. 

The luminosity and spectrum of radiation emitted during jet transition to sub-relativistic 
expansion are determined by the total energy E carried by relativistic plasma, by the num- 
ber density n of surrounding gas, by the fraction €b (e e ) of shock thermal energy carried 
by magnetic field (relativistic electrons), and by the shape of the electron distribution func- 
tion, which is commonly assumed to be a power law of index p = dlnn e /dln7 e m 2, 
where 7 e is the electron Lorentz factor (FWKOO). It is natural to assume that the values 
of €b, € e , and p are universal, since they are determined by the microphysics of the colli- 
sionless shock driven into the surrounding gas. Indeed, the constancy of p and e e among 
different bursts is strongly supported by observations. In bursts where p can be deter- 
mined accurately (e.g. Galama et al. 1998a, FWKOO, Stanek et al. 1999) p = 2.2 ± 0.1 
is inferred, a value consistent with numeric and analytic calculations of particle accelera- 
tion via the first order Fermi mechanism in relativistic shocks (Bednarz & Ostrowski 1998; 
Kirk et al. 2000; Achterberg et al. 2001). Universal values of p and e e are also inferred from 
the clustering of explosion energies (Frail et al. 2001) and of X-ray afterglow luminosity 2 
(Freedman & Waxman 2001; Berger et al. 2003a). The value of e B is less well constrained 
by observations. However, in cases where e# can be reliably constrained by multi waveband 
spectra, values close to equipartition are inferred (e.g. FWKOO). 



2 Apparently deviant values of p (Chevalier & Li 1999; Panaitescu & Kumar 2002) are inferred based on 
light curves, rather than spectra, and are sensitive to model assumptions (e.g. they depend on the assumed 
radial dependence of the ambient medium density). 
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The total energy E also appears to be universal. Although the apparent isotropic 7- 
ray energy, E 1; i so . varies by ~ 2.5 orders of magnitude between different bursts, a strong 
correlation between -E 7i i so . and 6j is observed, implying a narrow distribution of beam- 
ing corrected 7-ray emission, E 1 = d?E 1 j so j2 m 10 51 erg with roughly factor 3 spread 
(Frail et al. 2001). Moreover, the beaming corrected X-ray afterglow luminosity (at fixed 
time), Lx = #f£x,iso./2, which provides a robust estimate of the total kinetic energy car- 
ried by the relativistic plasma (Freedman & Waxman 2001), is also approximately constant 
(Berger et al. 2003a). Thus, variations in the radio flux during transition to sub- relativistic 
expansion are most likely due to variations in n, and possibly due to variations in e#. 

The prediction of L u ~ 10 3 V 10 QH z erg/sHz at ~ 1 yr delay was derived under the as- 
sumption of expansion into a uniform density medium with n ~ lcm -3 , consistent with the 
inferred values of n, which typically lie in the range of 10 5±1 cm~ 3 (e.g. 
(Bloom, Frail & Kulkarni 2003)). Similar predictions hold, however, for expansion into a 
wind with 

m = (M/lO^Moyr-^/^/lO^ms- 1 ) ~ 1 (Livio & Waxman 2000). rh ~ 1 is typi- 
cally adopted for modelling GRB afterglows in the scenario of expansion into wind (e.g. 
(Chevalier & Li 1999)), since the massive stars believed to be the progenitors of SNe Ib/c 
associated with GRBs (e.g. (Iwamoto et al. 1998; Woosley, Eastman & Schmidt 1999)), are 
observed to have winds with m > 1 (Willis 1991). For the case of SN1998bw, the density 
of gas surrounding the progenitor is constrained by the observed radio emission, which is 
consistent with synchrotron emission from a supernova shock wave propagating into a wind 
(Waxman & Loeb 1999; Li & Chevalier 1999). 

In order to address the question of whether or not an off-axis jet interpretation of 
GRB980425 is consistent with the long term radio observations of SN1998bw, we generalize 
in §2 the analysis of FWK00, which applies to expansion into a uniform medium, to the 
case of expansion into a wind. Our results are of general interest, beyond the analysis of 
GRB980425, since jet propagation into a wind, rather than into a homogeneous medium, 
may of course be characteristic for jets associated with SNe Ib/c in general. In §3 we discuss 
the implications of the results of §2 to the case of SN1998bw/GRB980425. Our conclusions 
regarding the nature of SN1998bw/GRB980425 and regarding the search for off-axis GRB 
signatures in nearby type Ib/c supernovae emission are summarized in §4. 

2. Model 

We first discuss in § 2.1 the hydrodynamics. Synchrotron emission is discussed in § 2.2. 
Some aspects of the model applicability are discussed in § 2.3. 
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2.1. Dynamics 

Let us consider a conical jet of opening angle 8j expanding into a r~ 2 density profile 
created by stellar mass loss. We assume a density profile 

where M is the mass loss rate and v w is the wind velocity. As long as T 3> l/Oj, the 
transverse size of causally connected regions within the flow, r/T(r), is smaller than the 
transverse size of the jet, Of, and the jet evolves therefore as if it were a conical section 
of a spherical relativistic blast wave. At this stage, the flow may be well described by the 
Blandford-McKee (1976) self-similar solutions, where 

\8tt Krc 2 J ' y } 

.Else, is the energy the flow would have carried had it been spherically symmetric. The true 
energy carried by the relativistic plasma is given by E = 6 2 Ei so J 2 (assuming a two-sided 
jet). At a radius r = Rg where T(r) drops to l/Oj, T(r = Rg) = l/Oj, the transverse size of 
causally connected regions exceeds Qf and the jet starts expanding sideways. At this stage 
the jet Lorentz factor drops exponentially with r (Rhoads 1999), and on a time-scale ~ Rg/c 
the flow approaches spherical symmetry. Using eq. 2 we have 

9 6]E lso . 9Ev w 18 E 51 
R e = = — = 1.5 x 10 cm, 3 

8tt Kc 2 Mc 2 rh K ' 

and the transition to spherical expansion takes place over a time scale 

tg^ Rg/C= 1.7— JT. (4) 

m 

Here, m = (M /lO^M^r" 1 ) /(^/lO^ms- 1 ). 

Eq. 3 implies that the mass enclosed within r < Rg is comparable to E/c 2 , and therefore 
that as the flow approaches spherical symmetry it also becomes sub-relativistic. For shock 
radii R > Rg the flow approaches therefore the self-similar Sedov-von Neumann- Taylor 
solutions describing expansion of a spherical strong shock wave into a r~ 2 density profile 
(e.g. Chapter XII of Zel'dovich & Raizer 2002). In these solutions, the shock radius is given 
by 

/ F \ 1/3 

R = m {K t2 ) ■ ^ 
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Here 7 is the adiabatic index of the gas, and £(7) is a dimensionless parameter of order 
unity. Exact determination of £ (7) requires a numerical solution of the flow variable profiles. 
However, a straight forward generalization of the Chernyi-Kompaneets approximation for the 
Sedov-von Neumann- Taylor solutions in a homogeneous medium (e.g. Chapter I of Zel'dovich 
& Raizer 2002) to the case of a r~ 2 density profile yields an analytic approximation for £, 



m = \ 



(7 + l) 2 (7-l) l V3 , fi , 
2tt(77-5) J ' {) 



For 7 = 5/3, appropriate for sub-relativistic flow, we have £ = 0.73. Denoting by -Rsnt the 
radius at which the shock velocity equals c, we have 

^ = ^£3^ = 0.38x10^ cm (7) 
9 S Mc 2 m K J 

We define £ S nt as the time at which the Sedov-von Neumann- Taylor solution gives R = c, 

E 

t S NT = 2/?snt/3c = 0.27— yr. (8) 

m 

Comparing Eqs. 3 and 7, we infer that the flow becomes sub-relativistic as it approaches 
spherical symmetry. 

The wind density at r = Rg is 

n e = 0.12—5- cm" 3 . (9) 

Eq. 9 shows that for rh ~ 1, the density n e is similar to that typical for the inter-stellar 
medium. Thus, the predicted radio signatures for a jet expanding into a rh ~ 1 wind are sim- 
ilar to those for a jet expanding into a typical inter-stellar medium (Livio & Waxman 2000). 
Rg in the case of expansion into uniform medium of density n = lnocm -3 is given by 
(Livio k Waxman 2000) 

*=(!i^T = ™ xwl \^T™> (io) 

and the corresponding time scales are (Livio & Waxman 2000) 

■ p \ i/3 / P \ 1/3 

-^51 \ ^ nnn ^51 



^ = 1.0 — yr, t SNT = 0.20 yr. (11) 

n J V n ' 
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2.2. Synchrotron emission 

At times t ^> £snt the flow is well described by the spherical non-relativistic self- 
similar solution, where the shock radius is given by Eq. 5. The long term afterglow observed 
in GRB970508 (FWKOO) and several other afterglows (Frail et al. 2003) can be explained 
as synchrotron emission from electrons accelerated by the collisionless shock to relativistic 
energy We assume that a constant fraction e e (e#) of the post-shock thermal energy is 
carried by relativistic electrons (magnetic field), and that the electron distribution function 
follows a power-law, dlnn e /dhi7 e = p for 7 e > j m . For simplicity, we assume that the 
shocked plasma is concentrated into a thin shell behind the shock, R/t] <C R, within which 
the plasma conditions are uniform (the Chernyi-Kompaneets approximation). This implies, 
in particular, rj — (7 + l)/(7 — 1). 

Under the above assumptions, the scaling of magnetic field amplitude B and of 7 m with 
time is given by (compare to Eq. A2 in FWKOO) 

b oc r\ 7m oc r 2 / 3 . (12) 

The flux and spectrum of synchrotron emission is then given by Eqs. A3 to A12 of ap- 
pendix A of FWKOO, with power-law indices of the t dependence in Eqs. { A6,A7,A8 }, 
{11/10, 1 — 3p/2, —3}, replaced with {11/6, —1 — 7p/6, —7/3} (due to the difference in tem- 
poral scalings between Eq. 12 and Eq. A2 of FWKOO). These Eqs. imply, in particular, that 
the characteristic synchrotron frequency of electrons with 7 e = / ~f m is given by 

v m = 2.3(3e e ) 2 (3e B ) 1 / 2 — [ — ) GHz, (13) 

and that the specific luminosity at frequencies v v m (and below the cooling frequency) is, 
for p — 2, 

• 9/4 —1/2 / + \ — 3/2 

L^l.Tx^W^^)- (JL) erg / S H, (14) 

Here too, the luminosity predicted for the case of expansion into rh — 1 wind is similar 
to that predicted for expansion into a uniform density medium (Livio & Waxman 2000), 

-1/2 / f \ - 9 / 10 

L v = 6.6 x 10 30 (3 ee )(36 B ) 3 / 4 nf (7^) (t^J er g/ sHz - ( 15 ) 
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2.3. Applicability and robustness 

The flow approaches spherical symmetry on a time scale t$ (Eq.4), and since £snt < tg 
(Eq. 8) it becomes (mildly) sub-relativistic at £snt < t <t$. Thus, an off-axis observer lying 
on a line of sight which makes a large angle, 9 o a ~ 1 rad, with the jet axis, will detect a 
flux comparable to that given by the model, Eq. 14, at a time £snt < to.a. < te, 

*o.a. ~l^yr. (16) 
m 

Observers located closer to the line of sight will detect a higher radio flux at earlier times, 
t < i .a.- At later times, t ^> tsNT, the flow approaches self-similarity and Eq. 14 provides a 
progressively more accurate approximation to the observed flux. On the other hand, the exact 
shape of the light curve at earlier time, t < t Q .a., is highly model dependent: The sideways 
expansion and the deceleration of the jet depend on the spatial distribution within the jet 
of the energy density and the Lorentz factor. These distributions are poorly constrained by 
current observations. Moreover, for given energy density and Lorentz factor distributions, 
an accurate calculation of jet expansion and deceleration can only be carried out numerically 
(e.g. (Ayal & Piran 2001; Granot et al. 2001)). 

Granot & Loeb (2003), for example, have recently considered emission from a point 
source on the jet axis observed by an off-axis observer in the case of expansion into a uniform 
density medium, with a simplified deceleration model. For E = 10 51 erg and n = lcm~ 3 
they find that an observer at o a . ~ 1 rad detects a flux similar to that detected by an 
on-axis observer at a time t _ a , = 0.4 yr. The specific luminosity they obtain at this time 
is 1.2 x 10 29 erg/sHz at 43 GHz, assuming {e e = 0.1, = 0.01, p = 2.5}, which (using the 
dependence of L v on e e , e# given by Eq. 14) corresponds to 8.0 x 10 30 erg/sHz at 10 GHz for 
e e = €b = 1/3. These results are similar to the analytic predictions of Eqs. 11, 15. 

In the Granot & Loeb (2003) analysis, the flux rises at t < t oa . nearly like a step 
function, and reaches at t somewhat larger than t Q . a . a maximum, which is an order of 
magnitude higher than the flux detected by an on-axis observer. This behavior is due (as 
also pointed out by the authors) to the simplified analysis, where only emission from a point 
source on the jet axis is considered. In a more realistic model, which includes emission from 
the entire jet and where sideways expansion and deceleration are calculated in more detail, 
the flux detected by a # o a . ~ 1 rad observer rises gradually at t < t Q . a . and the peak flux is 
expected to be (a few times) lower (e.g. (Granot et al. 2002)). 
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3. Implications to SN1998bw/GRB980425 

For fn ~ 1, as expected for GRB SN progenitors (see discussion following Eq. 9), an 
observer at large, 9 _ a , ~ 1 rad, offset is expected to measure at t Q . a . ~ f yr (Eq. 16) a specific 
luminosity L u ~ 10 30 erg/sHz at 10 GHz (Eq. 14, assuming e e ~ e B ~ 1/3). For a burst at 
cosmological distance this implies f u ~ 0.1 mJy This prediction is similar to that obtained 
for expansion into a homogeneous medium with density typical to that of the inter-stellar 
medium, n ~ O.lcm^ 3 (Eq. 15). Specific luminosities (fluxes) consistent with this prediction 
have been measured for GRB970508 (FWK00) and for several other long duration radio 
afterglows (Frail et al. 2003). 

Figure 1 compares the observed radio flux of SN 1998bw with the synchrotron flux 
predicted by the model described in § 2. The upper set of solid curves corresponds to 
fa = 1 (the peak in the 1 GHz light curve is due to self-absorption). The radio luminosity 
of SN1998bw/GRB980425 at 1 yr delay is 2.5 orders of magnitude lower than this model 
prediction. An off axis jet with "standard" energy E = 10 51 erg and "standard" e e ~ 1/3 
expanding into a wind with fa = 1 can therefore be ruled out under the assumption e B ~ 1/3. 
An off-axis jet interpretation of GRB980425 may be consistent with observations if it is 
assumed that either the wind of the progenitor of SN1998bw was atypical, with m < 1, or 
that the magnetic field energy density fraction has been atypically low, e B <C 1 (or both). 
Lower values of fa and oie B reduce the specific luminosity at the transition to sub-relativistic 
expansion, which is proportional to L v oc m 3//4 e^ //4 (Eqs. 14,8). A low value of m < 1 also 
delays the time at which the flow approaches sub-relativistic expansion (Eq. 16), and hence 
the time at which the radio emission approaches that of Eq. 14 for an off-axis, # . a . ~ 1 rad, 
observer, to ^> 1 yr. 

The radio emission of SN1998bw is consistent with synchrotron emission from electrons 
accelerated by the shock wave driven into the wind by the supernova ejecta (Waxman & Loeb 1999; 
Chevalier & Li 1999). Radio observations can therefore be used to constrain fa. As usually 
is the case for radio SNe, the data are not sufficient for determining all model parameters. In 
particular, there is a degeneracy in model predictions, that can be tested by observations, be- 
tween fa and e B , the fraction of the supernova post-shock thermal energy carried by magnetic 
field. Interestingly, the radio light curves imply that either fa or e B are atypically low: Data 
are consistent, e.g., with near equipartition value of e B and fa ~ 0.1 (Chevalier & Li 1999), 
and also with fa ~ 6 and e B ~ 10~ 6 (Waxman & Loeb 1999; Chevalier & Li 1999). 

If we adopt a model with e B <C 1 and fa > 1, we would need to assume e B < 10~ 4 1 
in order to reconcile radio observations with the predictions of an off-axis jet. In this case, 
the radio flux of the decelerated jet becomes undetectable. If, on the other hand, we adopt 
a model with l B near equipartition and fa ~ 0.1, the off-axis jet interpretation would be 
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consistent with observations for e# near equipartition. In this case, the transition to sub- 
relativistic expansion takes place on ~ 10 yr time scale, and assuming that we are at ~ 1 rad 
offset from the jet axis, the radio flux from the decelerated jet could not have been observed 
during the ~ 1 yr observations of SN1998bw. This is illustrated in Fig. 1 by the lower set of 
solid curves which describe model predictions for m — 0.1. 

Several clarification comments should be made at this point. The off-axis jet interpreta- 
tion of GRB980425 typically assumes that our line of sight is a few degrees away from the edge 
of a sharp-edged conical jet (Yamazaki, Yonetoku & Nakamura 2003; Granot et al. 2002). 
For a jet with sharp edges, the flux drops rapidly when the line of sight makes an angle with 
the jet axis which is larger by ~ 1/T ~ 0.01 then the jet opening angle, which constrains 
our line of sight to be within a few degrees of the jet edge. This interpretation would be 
inconsistent with observations for the case of e B ~ 1/3 and m ~ 0.1, since in this case jet 
deceleration would bring our line of sight into the jet's radiation "beaming cone", leading to 
strong radio emission on time scale < 1 yr. Our line of sight must make a large angle with 
the jet axis in order to avoid observing the strong radio emission from the decelerated jet at 
t < 1 yr. The observed 7-ray flux of GRB980425 can be explained in this case by assuming 
that the jet is not sharp-edged, but rather has "wings" that extend to ~ 1 rad and produce 
the observed low 7-ray luminosity. Alternatively, Compton scattering of photons into our 
line of sight may allow a large off- axis orientation (Eichler & Levinson 1999). 

For m = 0.1, the time at which an off- axis, oa . ~ 1 rad, observer detects a flux 
comparable to that predicted by Eq. 14 is t Q . a . ~ 10 yr (Eq. 16). At this delay, the predicted 
radio flux is a few mJy, and the predicted optical flux is a few /iJy (my = 22.5), assuming 
e B ~ 1/3. What is the current, t ~ 5 yr, flux predicted by this model? As discussed in 
§ 2.3, the detailed time dependence of flux observed at t < t Q . a . is highly model dependent. 
In realistic jet models, however, the flux is not expected to vary strongly at t ~ t Q . a .. Hence, 
a ~ 1 mJy 10 GHz flux and m v ~ 23 optical flux are expected in this model. 

4. Discussion 

Simple analytic expressions are given in § 2 for the specific luminosity emitted by a 
GRB fireball expanding into a wind, after it had decelerated to sub-relativistic speed. At 
this stage radiation is emitted roughly isotropically. An off-axis observer lying on a line 
of sight which makes a large angle, # oa . ~ 1 rad, with the jet axis, is predicted to detect 
a flux comparable to that given by the model, Eq. 14, at a time t ~ t Q .a. = lEmm^ 1 yr, 
where E = 10 51 i?5i erg. The specific luminosity at 10 GHz at t ~ t Q _ a , is L v ps 0.24 x 
10 30 (3e e )(3es) 3 / 4 m 3 / 4 £ , 5ierg/s Hz. At later times the model provides a progressively more 
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accurate approximation to the observed flux. The exact shape of the light curves at earlier 
time, t < t o a ., is highly model dependent. In particular, it depends on the (unknown) spatial 
jet structure. 

The low radio luminosity of SN1998bw, compared to that expected from a decelerated 
GRB jet at ~ 1 yr delay, may be consistent with the off-axis jet interpretation of GRB980425 
provided that either the magnetic field energy fraction is atypically low, e# < 10~ 4 , or 
that the density of the wind surrounding the progenitor is lower than typically expected, 
fn = (M/lO~ 5 M yr^ 1 )/(f w /lO 3 kms _1 ) ~ 0.1. Lower values of fn and of e# reduce the 
specific luminosity at the transition to sub-relativistic expansion. A low value of m <C 1 
further delays the time at which the flow approaches sub-relativistic expansion, as illustrated 
in Fig. 1. In the former scenario, e B < 1CT 4 , the flux from the decelerated jet becomes 
undetectable. We consider this scenario less likely, however, since we expect e^, which 
is determined by the shock micro-physics, to be similar for different bursts, and es close 
to equipartition is inferred from radio observations of other bursts. The latter scenario, 
fn ~ 0.1, is consistent with the constraints imposed on fn by the observed radio emission 
from the supernova shock driven into the wind. In this scenario, transition to sub-relativistic 
expansion is expected over ~ 10 yr time scale. A ~ 1 mJy 10 GHz flux and my ~ 23 optical 
flux are expected in this model at present, and the angular source size is expected to be 
~ 10 mas. 

Our analysis also demonstrates that in order to place robust constraints on the fraction 
of SNe Ib/c associated with cosmological GRBs, radio follow up of such SNe may need to 
be carried over a multi-year period, since for lower values of fn strong radio emission may 
be detected by an off-axis observer only after several years (Additional motivation for multi- 
year monitoring is the possibility to identify radio emission from "failed GRBs," with low 
velocity and large baryon load jets (Totani 2003)). The values of n inferred from afterglow 
observations typically lie in the range of 10 a5±1 cm~ 3 (e.g. (Bloom, Frail & Kulkarni 2003)). 
As explained in § 2.1 and § 2.2 (see also (Livio & Waxman 2000)), the predicted long-term 
radio signatures for a jet expanding into am~l wind are similar to those for a jet expanding 
into a typical inter-stellar medium n ~ lcm -3 . For n ~ lcm~ 3 , or fn ~ 1, strong radio 
emission should indeed be observed by an off axis observer at t ~ 1 yr (see Eqs. 16, 11). 
However, observations of GRB980425/SN1988bw demonstrate that, at least in some fraction 
of the cases, strong radio emission from a decelerated jet would be detected by an off-axis 
observer only over a longer time scale. 

Finally, it should be kept in mind that the wind mass-loss rate from the progenitor of 
SN1998bw may have been time dependent. The wind density profile at distances < 1 pc 
corresponds to the wind history over ~ 10 3 yr preceding the explosion (for a wind speed 
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of 10 3 km/s). We have no direct information on the steadiness of the mass loss rate from 
massive stars so close to the end of their evolution. A higher mass loss rate at earlier times, 
closer to the typically expected m ~ 1, would lead to a higher radio flux than predicted in 
figure 1 at late times. 

I thank B. Paczyhski for discussions that initiated the study described in this manuscript. 
EW is partially supported by AEC grant and a Minerva grant, and is an incumbent of the 
Beracha Career development Chair. 
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Fig. 1. — The radio flux of SN 1998bw compared with predictions of the model described 
in § 2 for a jet of energy E = 10 51 erg expanding into a wind. The upper, lower sets of 
solid curves show model flux at t > tgNT for m = 1, m = 0.1 respectively. Equipartition, 
e e = ts = 1/3, and p = 2.2 are assumed. An off-axis observer lying on a line of sight which 
makes a large angle, o a . ~ 1 rad, with the jet axis, is predicted to detect a flux comparable 
to that given by the model, Eq. 14, at a time t ~ i . a . = \mT x yr. At later times the model 
provides a progressively more accurate approximation to the observed flux. On the other 
hand, the exact shape of the light curves at earlier time, t < i . a ., is highly model dependent. 
In particular, it depends on the (unknown) spatial jet structure. The dashed curves, added 
for illustration, follow the temporal dependence found in the numerical simulations described 
in Granot et al. (2002). 



